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We investigate heavy-ion multinucleon transfer reactions using the coupled-channels formalism. 

We first use the semi-classical approximation and show that a direct coupling between the entrance 
and the pair transfer channels improves a fit to the experimental one- and two-neutron transfer cross 
sections for the ^®Ca-|-®®Zr and ®°Ni-|-^^®Sn systems. We then discuss the validity of the perturbative 
approach and highlight the effect of high order terms. The effect of absorption is also investigated 
for energies around the Coulomb barrier. Finally, we use a quantal coupled-channels approach to 
achieve a simultaneous description of the fusion cross sections and the transfer probabilities for the 
^°Ca-|-®®Zr reaction. We find a significant effect of the couplings to the collective excited states on 
the transfer probabilities around the Coulomb barrier. 
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I. INTRODUCTION 

A transfer of a few nucleons during the reaction pro¬ 
cess plays an important role in heavy-ion sub-barrier fu¬ 
sion reactions [T]. The coupling to a transfer degree of 
freedom is in general weaker than couplings to collective 
excitations. However, when the coupling is sufficiently 
large, e.g., for a transfer channel with a positive Q'Value, 
it is expected that the transfer coupling leads to an extra 
enhancement of fusion cross sections below the Coulomb 
barrier HIE]. 

It has been known well that two-neutron transfer reac¬ 
tions provide a unique tool to study the pair correlation 
between nucleons [I9U22] . Recently two experiments for 
heavy-ion transfer reactions have been carried out at en¬ 
ergies far below the Coulomb barrier [231 [21] in order to 
study the relation between the transfer of one neutron 
and that of a pair of neutrons. See also Refs. [JSl l26] . 
The experimental data show an enhancement of the two- 
neutron transfer cross sections compared to a simple esti¬ 
mate based on the independent picture, that is, a square 
of the transfer probability of one neutron. 

Several theoretical studies have been performed in or¬ 
der to understand the reaction dynamics of the pair 
transfer process. Those include the time-dependent per¬ 
turbation theory based on the semi-classical approxi¬ 
mation m uni uzHiH], the second-order distorted wave 
Born approximation (DWBA) [22l [301 - 132] , and the semi- 
classical [26] and the quantal mis] coupled-channels ap¬ 
proaches. Recently, the time-dependent density func¬ 
tional theory (TDDFT) has also been employed to inves¬ 
tigate the transfer reaction. An advantage of this method 
is that the formalism takes into account simultaneously 
the structure and the dynamics. Although this method 
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reproduces the individual transfer of nucleons |33H36| . 
the enhancement factor of a two-neutron transfer prob¬ 
ability is underestimated with this method |37j even if 
the pairing correlation is taken into account [38j . This 
implies a necessity to go beyond the mean-field dynamics 
with pairing in order to take into account all the remain¬ 
ing correlations |39H42]. 

In this paper, we carry out a phenomenological study 
in order to understand the transfer dynamics using the 
formalism of the semi-classic and the quantal coupled- 
channels approaches. Our aim in this paper is twofold. 
One is to understand the reaction dynamics of the two- 
neutron transfer reactions. We emphasize that the reac¬ 
tion dynamics is so complicated that many microscopic 
approaches still lack of prediction for the two-neutron 
transfer. In particular, the nature of the two-neutron 
transfer reaction, that is, the relative importance between 
the direct and sequential processes, is still under discus¬ 
sions 1311211 HZl 1^ ISi 133H1S] . The second aim of this 

paper is to describe subbarrier fusion reactions using the 
transfer coupling form factors which are consistent with 
the transfer cross sections. For these purposes, we phe¬ 
nomenologically adjust the parameters in the form fac¬ 
tors for the transfer couplings, rather than computing 
them microscopically. 


The paper is organized as follows. In Sec. II, we em¬ 
ploy the semi-classical time-dependent coupled-channels 
method and discuss the sequential and the direct natures 
of the two-neutron transfer process. In Sec. Ill, we use 
the transfer coupling form factors obtained in Sec. II to 
discuss the role of absorption in the two-neutron trans¬ 
fer reactions. In particular, we investigate the interplay 
between the transfer and the fusion processes using the 
quantal coupled-channels approach. We then summarize 
the paper in Sec. IV. 
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II. SEMI-CLASSICAL APPROACH 


''°Ca+s®Zr 


“Ni+ii®Sn 


A. Transfer probability 


Since the experimental transfer cross sections are often 
analyzed using the semi-classical method [TS], we first 
employ it to investigate the nature of the two-neutron 
transfer. 

With the semi-classical approximation to coupled- 
channels equations, one assumes a Rutherford trajec¬ 
tory for the relative motion r{t) between the colliding 
nuclei. This yields a time-dependent field for the in¬ 
trinsic motions in the projectile and the target nuclei. 
That is, the nuclear intrinsic wave function is expanded 
as l^'(t)) = Cn{t)\n), n = 0 corresponding to the 

entrance channel, and it is then evolved in time as 

( 1 ) 

n' 

with Cn{—oo) = Sn,o- Here, Hnn’it) is given as 

n 

nn'(0 — ^n^n,n' ^nn' [^W], (2) 


where is the excitation energy for the channel n and 
Vnn'[r(t)] is the coupling form factor evaluated along the 
classical trajectory, r{t). The probability for the channel 
n is given as = |c„p at t = -|-c». 

For the transfer problem, n corresponds to the num¬ 
ber of transferred nucleons. By truncating the transfer 
channels at u = 2, the coupling Hamiltonian T-L reads, 

/ 0 Voi{t) Vo2{t) \ 

n{t) = Volit) -Qi Hi 2 (t) , (3) 

\ Vo2{t) Vi2{t) —Q 2 / 


where e„ = —Qn is the transfer Q-value for each parti¬ 
tion, n. For the coupling form factor Vnn', we employ 


Vnn' (r) 



1^1 _|_ 


(4) 


This function has the exponential form. 


V„n' (r) 


•\/47r 


g-(r-rp)/a„ 


(5) 


for r ^ Tp. Note that this parametrization differs from a 
frequently used form, 14 „'(r) = where /(r) 

is a Fermi function, only for small values of r. We find 
that the solutions for the coupled-channels equations are 
numerically more stable with the parametrization given 
by Eq. Q. 

In the actual calculations presented below, we use 
Tp = 1.1 X fni; where At and Ap are the 

mass numbers for the target and the projectile nuclei, 
respectively. We start the time-dependent Schrbdinger 
equation, Eq. Q, at at an initial distance nni = 30 fm 
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FIG. 1: The ground state energies for the one-nentron (In) 
and the two-neutron (2n) transfer processes for the ^°Ca-|-®®Zr 
and ®°Ni-|-^^®Sn systems. The dashed lines denote the opti¬ 
mum Q-value for each transfer channel. 


where the coupling is negligible. Following Ref. [29], 
we evaluate the time dependence of the coupling form 
factor Vnn'{t) by averaging the two Rutherford trajec¬ 
tories, r„(t) and as r{t) = (r„(t) -I- r„/(t)V2, 

where r„(t) and are evaluated with the energy of 

En = Pcm + Qn and p„/ = Ecm + Qn' , respectively. In 
order to avoid a dependence of the results on the initial 
position, Tini, all the trajectories are arranged so that 
the distance of the closest approach is reached at the 
same time, tmin- We follow the time evolution until the 
two fragments are separated with a distance of 30 fm. 
We have checked that the coupling among the channels 
are practically negligible after this distance, and thus the 
transfer probabilities are not altered. 

In order to reduce the number of adjustable parame¬ 
ters, we employ two different schemes. The first scheme 
corresponds to a pure sequential process, for which we set 
V 02 to be zero while Vbi and F 12 are allowed to be differ¬ 
ent. On the other hand, the second scheme corresponds 
to a direct two-neutron transfer process, for which we set 
Voi{r) = Vi 2 {r) but Vb 2 7 ^ 0 in general. The enhance¬ 
ment of two-neutron transfer process is explained in a 
different way with these two schemes. In the pure se¬ 
quential transfer scheme, the enhancement factor is due 
to the increase of the probability for the transfer of the 
second neutron after the hrst neutron is transferred, that 
iS) IH 12 I > |Voi|- On the other hand, in the direct transfer 
scheme, the enhancement is due to the additional cou¬ 
pling between the entrance and the two-neutron transfer 
channels. 

We test these schemes on the two recent experiments 
for the 40Ca-f®®Zr and 60Ni-7“6Sn systems |2ll|24|- In 
these references, the experimental data, taken at back¬ 
ward angles in the center of mass frame, are given in 
terms of the transfer probability dehned by the ratio of 
the transfer cross sections to the Rutherford cross sec¬ 
tions, = da^n/dan, with a; = 1, 2 , • • •, as a function 
of the distance of the closest approach, D, for the Ruther¬ 
ford trajectory. Figure summarizes the ground state 
energies for these systems. The Q-value for the ground 
state to the ground state transition, Qgg, is positive both 
for the one-neutron and the two-neutron transfer chan- 
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FIG. 2: (Color online) The transfer probabilities for the 
"'^Ca+^^Zr (the upper panel) and ®“Ni+^^®Sn (the lower 
panel) reactions as a function of the distance of the closest 
approach, D, for the Rutherford trajectory. The experimen¬ 
tal data, taken form Refs. | 231124| . for the one-neutron (the 
crosses) and the two-neutron (the squares) transfers are com¬ 
pared to transfer probabilities for the ground state to the 
ground state transitions obtained with the pure sequential 
scheme (the thin solid line for the In transfer and the dotted 
line for the 2n transfer) and with the direct transfer scheme 
(the thick solid line for the In transfer and the dashed line 
for the 2n transfer). 

nels for the "^°Ca-|-®®Zr system, while it is negative for 
the one-neutron transfer and positive for the two-neutron 
transfer for the ®°Ni-|-^^®Sn system. 

Let us first consider a case with the ground state to the 
ground state Q-value, that is, a case with Qn = Qgg(n) 
in Eq. (^. We adjust the coefficients a„„' and /3nn' in 
Eq. @ by fitting the experimental data for D > 13.5 fm, 
for which the coupling is weak and only the first order 
dynamics is important. Figure shows a comparison of 
the calculated transfer probabilities so obtained with the 
experimental data. Those are obtained by varying the 
center of mass energy with fixed values of the scatter¬ 
ing angle, that is, dc.m. = 140 degrees for the ^°Ca-|-®®Zr 
and ®°Ni-|-^^®Sn systems. For each energy, the impact 
parameter is determined so that the scattering angle for 
the Rutherford trajectory for the entrance channel is con¬ 
sistent with the given value of the scattering angle. In 
the figure, the thin solid and the dotted lines denote the 
results with the pure sequential scheme for the In and 2n 
transfer processes, respectively. On the other hand, the 
thick solid and the dashed lines denote those with the 
direct transfer scheme. 

The values of the parameters are summarized in Ta- 


TABLE I: The parameters for the transfer coupling form fac¬ 
tors for the ^®Ca-|-®®Zr reaction. For all the rows, a\ and /3i 
are those for the coupling between the entrance and the In 
transfer channels. For the sequential scheme (seq.), 02 and /32 
are those for the coupling between the In and the 2n trans¬ 
fer channels, that is, 02 = ai 2 and d 2 = 012 ■ For the direct 
transfer scheme (dir.), 02 and 02 are those for the coupling 
between the On and the 2n transfer channels, that is, 02 = ao 2 
and 02 = 002- These parameters are obtained either with the 
ground state to the ground state transfer Q-values (Qgg) or 
with the optimum Q-values (Qopt). 


Q 

Scheme 

ai (fm) 

(MeVfm) 

02 (fm) 

02 (MeVfm) 

Qgg 

seq. 

1.015 

121.8 

0.726 

4470 

Qgg 

dir. 

1.309 

39.05 

0.727 

443 

Qopt 

seq. 

1.113 

80.33 

1.458 

107.4 

Qopt 

dir. 

1.230 

51.56 

0.700 

278.9 


TABLE II: Same as Table|^ but for the ®®Ni-|-^^®Sn reaction. 


Q 

Scheme 

oi (fm) 

(MeVfm) 

02 (fm) 

02 (MeVfm) 

Qgg 

seq. 

0.877 

138.9 

1.36 

131.9 

Qgg 

dir. 

1.18 

40.6 

0.602 

363.4 

Qopt 

seq. 

0.893 

126.6 

1.38 

109.4 

Qopt 

dir. 

1.13 

47.1 

0.600 

384.2 


bles I and II for the "^^Ca-l-^^Zr and ®°Ni-|-^^®Sn systems, 
respectively. Notice that the direct two-neutron transfer 
has a shorter range with small 02 than the sequential cou¬ 
pling, in accordance with the argument in Ref. |46j . One 
can also notice that with the sequential scheme one has 
to take a considerably different parameters for V 12 from 
those for Vqi in order to reproduce the experimental two- 
neutron transfer probability. This may appear unnatural 
since a drastic change of structure is not expected even 
after one neutron is transferred, although it may simply 
mock up a constructive interference of different transfer 
paths [17]. 

In Fig. an unexpected behavior can be seen for 
the "^°Ca-|-®^r reaction with the pure sequential scheme 
(the dotted and the thin solid lines). That is, the transfer 
probabilities oscillate as a function of I?, even though the 
exponential tail is correctly reproduced. This behavior 
is due to the large Q-value for the two-neutron transfer 
channel, (32=5.5 MeV. 

In reality, however, the transfer takes place mainly to 
excited states, rather than to the ground state [HI HHj. 
In fact, the optimum (3-value for a neutron transfer is 
Qopt = 0 [49|, and the coupling to the ground state is 
much weaker [^. Figure]^ is obtained by using the op¬ 
timum (3-values, that is, (3=min(0,(3gg) (see Fig. I^. 
See Tables I and II for the parameters. As expected, the 
agreement with the experimental data is much improved. 

In Fig. one can notice that the one-neutron trans¬ 
fer probability is well reproduced with the direct transfer 
scheme. In contrast, a significant deviation is seen for 
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FIG. 3: (Color online) Same as Fig. but with the optimum 
Q-values (see the text). 


small values of D with the sequential scheme. This can 
be attributed to the sequential nature of the transfer dy¬ 
namics. For large values of D, the transfer probabilities 
are small and the sequential transfer to the 2n channel 
does not modify Pi. However, as D decreases, the two- 
step process becomes significant, reducing the probabil¬ 
ity for the one-neutron transfer channel. With the direct 
two-neutron transfer scheme, on the other hand, the two- 
neutron transfer takes place mainly from the entrance 
channel (see Sec. II-C below), as the coupling between 
the In and the 2n channel is kept to be relatively weak. 

In order to better understand this phenomena, one can 
use the time-dependent perturbation theory by separat¬ 
ing the contribution of each order. Figure|^shows the ap¬ 
proximate solutions for the one-neutron transfer reaction 
for the ^°Ca-|-®®Zr system obtained with the sequential 
two-neutron scheme. Here, the coupling parameters are 
readjusted in order to reproduce the experimental data 
for large D with the first order dynamics. Up to the first 
order, the one-neutron transfer probability is found to be 
exponential with the distance of the closest approach in 
the whole range of D shown in the figure. As a conse¬ 
quence, the calculation appears to be consistent with the 
experimental data. This remains to be the same up to the 
second order, since there is no second order contribution 
to the one-neutron transfer channel. On the other hand, 
if one considers up to the third order, the result is dras¬ 
tically changed and becomes close to the full order result 
obtained by solving the time-dependent coupled-channels 
equations, Eq. without using the perturbation the¬ 
ory. That is, the third order process now includes the 
diminution of the one-neutron transfer probability due 
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FIG. 4: (Color online) The one-neutron transfer probability 
as a function of the distance of the closest approach for the 
"^°Ca-|-®®Zr reaction. The sequential two-neutron scheme with 
the optimum Q-value is employed. The dashed line, the trian¬ 
gles, and the dotted line are the results of the time-dependent 
perturbation theory with the first, the second, and the third 
orders, respectively. Those results are compared to the exact 
solution (the solid line) and to the experimental data (the 
crosses). The coupling parameters used here are ai =1.55 fm. 
Pi = 26.4 MeV fm, 02 =1.286 fm and P 2 =150.4 MeV fm. 


to the second order process to the 2n channel. From this 
comparison, we can conclude that the perturbative cal¬ 
culations must be used with caution for an application 
to multinucleon transfer calculations, since a good agree¬ 
ment with experimental data may be an artifact of the 
first order perturbation theory. 

Notice that the sequential scheme largely underesti¬ 
mates the transfer probabilities at small D. In this 
region, the absorption effect would be important, but 
the effect of absorption will always decrease the transfer 
probabilities, as we will discuss in Sec. H-C. Evidently, 
a pure sequential two-neutron transfer is not compati¬ 
ble with the experimental data, at least with this simple 
model. See also Refs. [51131] for a similar conclusion. We 
have confirmed that this conclusion remain the same even 
if we include a few states around the optimum Q-value 
for each transfer partition. 

Notice that the semi-classical transfer calculations 
have indicated that the two-neutron transfer process oc¬ 
curs predominantly with a sequential process 
The present result is somewhat in contradiction with the 
previous finding of the semi-classical method. The dif¬ 
ference may be due to the fact that the previous semi- 
classical calculations take into account many intermedi¬ 
ate In transfer channels with a wide distribution of exci¬ 
tation energy while we include only a single (or at most 
a few) channel. In this sense, our coupling scheme may 
have to be regarded as an effective one, which implicitly 
takes into account the effect of many other intermediate 
channels. It would be an interesting future work to clarify 
how the elimination of the intermediate states gradually 
changes the nature of the two-neutron transfer couplings. 
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FIG. 5: (Color online) The two-neutron transfer probability 
for the ^®Ca-|-®®Zr system as a function of the distance of the 
closest approach, D. The direct two-neutron transfer scheme 
with the optimum Q-values is employed. The results of the 
time-dependent perturbation theory up to the first and the 
second orders are shown by the dotted and the solid lines, re¬ 
spectively. The result of the second order calculation neglect¬ 
ing the direct two-neutron transfer couplings is also shown by 
the dashed line (see text). The experimental data are taken 
from Ref. | 23| . 


B. Nature of two-neutron transfer 


In the previous subsection, we have argued that the di¬ 
rect two-neutron transfer scheme is more plausible than 
the pure sequential transfer scheme as long as a simple 
coupling scheme as in Eq. ^ is employed. Let us then 
investigate the nature of the two-neutron transfer process 
by assuming the direct transfer scheme with the optimum 
Q-values. Figure shows the result for the probability 
of the two-neutron channel for the “^^Ca-l-^^Zr system ob¬ 
tained by the perturbation theory. To this end, we use 
the same parameters as those used in Fig. The dot¬ 
ted line denotes the result of the first order perturbation 
theory, which includes only the direct population of the 
2n transfer channel from the entrance channel. The solid 
line, on the other hand, denotes the result of the second 
order perturbation theory, which in addition includes the 
sequential two-neutron transfer via the intermediate In 
transfer channel. For a comparison, the figure also shows 
the calculation with the second order contribution only 
(the dashed line), that is, the calculation neglecting the 
direct two-neutron coupling. One can see that the se¬ 
quential process is negligibly small as compared to the 
first order process. That is, the two-neutron transfer in 
this model is largely dominated by the direct two-neutron 
transfer from the entrance channel partition. Notice that 
this conclusion is a direct consequence of the simplified 
transfer model with the choice of Vi 2 {r) = Voi(r), and it 
would need to be confirmed with a more general model. 


TABLE III: The parameters of the absorbing potential for 
each channel required to fit the experimental data. 


System 

Channel 

Wo (MeV) Rw (fm) 

aw (fm) 


On and In 

9.96 

11.88 

0.12 


2n 

9.89 

12.10 

0.18 

“Ni-k“®Sn 

On and In 

10.03 

12.01 

0.08 


2n 

10.08 

12.25 

0.16 


C. Effect of absorption 

While the experimental data for the two-neutron trans¬ 
fer are well described with the direct two-neutron transfer 
scheme for the distance of the closest approach of 13 fm 
or larger, a discrepancy is found for U < 13 fm (except 
for the ^°Ca-|-®®Zr system with Qgg). While the calcu¬ 
lations predict a nearly exponential behavior of D, the 
experimental data show a much slower dependence in this 
region. This could be attributed to several mechanisms 
beyond the present model, such as a deviation of the 
classical trajectory from the Rutherford trajectory, a loss 
of one-to-one correspondence between the experimental 
and the theoretical definitions of the transfer probabili¬ 
ties (that is, the quantity dUxn/d^^R may not be able to 
be interpreted as a transfer probability at small values 
of 13), and a deviation of the coupling form factor from 
that given by Eq. (|^. 

In the following, we choose to test a hypothesis that 
this effect originates from the absorption of the wave 
function, which corresponds to a capture and/or inelas¬ 
tic excitations. In order to take into account this effect 
in a simple way, we add an imaginary potential to the 
diagonal part of the Hamiltonian, Eq. (|^, with the ex¬ 
pression. 


iW{r) 


-iWo 

1 -I- exp[(r — Rw)/o.w] ’ 


( 6 ) 


for which the parameters Hq, Rwj and aw may be dif¬ 
ferent for each channel. We adjust these parameters in 
order to reproduce the experimental transfer probabili¬ 
ties using the direct two-neutron transfer scheme with 
the optimum Q-values. The resultant values for the pa¬ 
rameters are summarized in Table III, and the fit to the 
experimental transfer probabilities is shown in Fig. 
One can see that a good agreement between the calcula¬ 
tions and the experimental data is achieved, including the 
region with small values of the distance of the closest ap¬ 
proach, although the good reproduction of the data may 
be due to the large number of adjustable parameters. 

The single-channel calculations with those absorbing 
potentials lead to the survival probability for each chan¬ 
nel shown in Fig. Notice that this is not the sur¬ 
vival probability for the whole transfer process, that is, 
Psurv = I Cop + |cip -I- |c 2 p, but the single-channel prob- 
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FIG. 6: (Color online) The transfer probabilities as a func¬ 
tion of the distance of the closest approach, D, calculated by 
including the effect of absorption of the classical trajectories 
with the imaginary potential. The direct two-neutron trans¬ 
fer scheme with the optimum Q-value is employed. The solid 
lines denote the results for the one-neutron transfer process, 
while the dashed lines are for the two-neutron transfer pro¬ 
cess. The experimental data are taken from Refs. |23l I24| . 


TABLE IV: The parameters for the complementary error 
function given by Eq. Q which fit the results of the time- 
dependent calculations for the survival probability for each 
channel. 


System 

Ghannel 

B (MeV) 

a (MeV) 

^°Ca-b®®Zr 

On 

95 

1.3 


In 

95 

1.3 


2n 

92 

1.8 

60Ni-^ii6sn 

On 

166 

1.3 


In 

166 

1.5 


2n 

160 

2.8 


ability corresponding to that given by [H] 

Psnrl\E) = 1 - exp Q ^ W[r^{t)]dt^ . (7) 

We find that the survival probability for each channel, 
Psu ^, can be well parametrized with the complementary 
error function, 

«(i5) ~ i + ierfc (^) , (8) 

with the parameters given in Table IV (see Fig. [^. No¬ 
tice that the values of the parameters are considerably 
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FIG. 7: (Color online) The s-wave survival probability for 
each channel as a function of the center of mass energy. It is 
obtained by solving the time-dependent Schrodinger equation 
with a single channel including the imaginary potential. The 
crosses and the squares show a fit with the the complementary 
error function given by Eq. ([^. 

different between the two-neutron transfer channel and 
the one-neutron transfer channel. The meaning of the 
imaginary potential is to take into account the loss of 
flux from the model space explicitly taken into account 
in the calculations to the outside, which include the tun¬ 
neling through the barrier (that is, the capture) and the 
inelastic processes. In either case, the imaginary poten¬ 
tial reduces the transfer probabilities. The different val¬ 
ues of the parameters for each channel implies that the 
effect of absorption may be different in each channel. We 
find that it is essential to have the channel dependence 
in the parameters of the absorbing potential in order to 
account for the diminution of the 2n transfer probability 
without significantly degrading the reproduction of the 
In transfer channel. 


III. QUANTAL COUPLED-CHANNELS 
APPROACH 

We have shown in the previous section that the role 
of absorption may differs among the different transfer 
channels. We have achieved this conclusion using the 
semi-classical method. A drawback of the semi-classical 
method is that the tunneling process is not easy to be 
described with it, although one may still do it using the 
time-evolution along the imaginary time axis [Sni[5T]. For 
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this reason, in this section, we shall use the full quantal 
coupled-channels approach, which can also be applied to 
a simultaneous description of fusion cross sections and 


transfer probabilities. 

The coupled-channels equations for a total angular mo¬ 
mentum J read 1501. 


2^ dr^ 


ZpZre^ 

„ ’ + VNir) + iWN{r) + - 

r 


-E 


Un{r) + ^ Vnm{r)u^{r) = 0 , 


(9) 


TABLE V: The parameters for the nuclear potential used in 
the coupled-channels calculations for each system. 


System 

Vo 

ro 

flO 

Wo 

rw 

aw 


(MeV) 

(fm) 

(fm) 

(MeV) 

(fm) 

(fm) 


87.0 

1.13 

0.700 

40 

1.1 

0.2 

®°Ni+“®Sn 

80.0 

1.1 

0.487 

20 

0.9 

0.487 


where /i is the reduce mass, E is the center of mass en¬ 
ergy, and Zp and Zp are the charge numbers of the pro¬ 
jectile and the target, respectively. We have used the 
iso-centrifugal approximation [50] and assumed that the 
angular momentum does not change for each channel. 
We use a Woods-Saxon parametrization for the nuclear 
potential, Vm, that is. 


VN{r) 


-Vo 

1 -hexp[(r- i?o)/ao] 


( 10 ) 
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with i?o = ro(Ap^ -|- We use the Woods-Saxon 

parametrization also for the imaginary potential, Wjv, as 
in Eq. (|^. 

We solve the coupled-channels equations using a mod¬ 
ified version of the computer program CCFULL [52] in or¬ 
der to construct the transfer cross sections from the S- 
matrix. The transfer probability is then calculated using 
the same definition as the experimental probability, that 
is, = da xn I dap. 


A. Quantum effect 

Let us first solve the coupled-channels equations in¬ 
cluding only the transfer channels, and examine the va¬ 
lidity of the semi-classical trajectory approximation. To 
this end, we use the direct two-neutron transfer scheme 
with the optimum Q-values, for which the parameters for 
the transfer coupling form factors are given in Tables |T| 
an m We use the parameters in Table [Vj for the nucleus- 
nucleus potential (both for the real and the imaginary 
parts). In order to simplify the comparison between the 
semi-classical and the quantal calculations, we neglect 
the effect of absorption on the transfer probabilities. For 


FIG. 8: (Color online) A comparison between the trans¬ 
fer probabilities obtained with the semi-classical approxima¬ 
tion (the crosses and the circles) and those with the quantal 
coupled-channels calculations (the solid and the dashed lines) 
for the "‘°Ca-|-®®Zr and ®°Ni-|-^^®Sn systems. 


this purpose, the real and the imaginary parts of the 
nuclear potential are not taken into account in the semi- 
classical calculation. Moreover, we choose the parame¬ 
ters of the nuclear potential so that they yield a higher 
barrier than the systematics. 

Figure [^ shows a comparison between the quan¬ 
tal coupled-channels calculations and the semi-classical 
coupled-channels calculations for the same coupling po¬ 
tentials. One can find a good agreement between the 
two calculations for the distance of the closest approach 
oi D > 13.5 fm, that is, the two calculations differ only by 
about 10 %. The difference between the two calculations 
is more significant for D < 13 fm, where the quantum ef¬ 
fect as well as the absorption process play an important 
role. 

Evidently, the conclusions obtained in the previous sec¬ 
tion with the semi-classical method remain unchanged 






















even if we use the quantal coupled-channels method. 


B. Role of absorption 



12.5 13.0 13.5 14.0 14.5 15.0 


D [fm] 

FIG. 9: (Color online) The one-neutron transfer (the upper 
panel) and the two-neutron transfer (the lower panel) proba¬ 
bilities as a function of the distance of the closest approach for 
the ®°Ni-|-^^®Sn reaction. The solid and the dotted lines are 
obtained with the nuclear potential which yields the barrier 
height of B=166 MeV and 179.9 MeV, respectively. 

We next discuss the role of absorption in the transfer 
reactions by fully taking into account the tunneling effect. 
For this purpose, we still use the three-channel problem 
in the previous subsection. Figure [^compares two calcu¬ 
lations for the transfer probabilities for the ®°Ni-|-^^®Sn 
system. The first calculation shown by the dotted line 
is the same as the dashed line in Fig. This result is 
obtained with the nuclear potential of Table V, which 
yields the barrier height of i? = 179.9 MeV. Notice that 
the Broglia-Winther potential [55] leads to a somewhat 
lower barrier height, that is, B = 166 MeV. We therefore 
repeat the transfer calculations using the potential with 
Vo=116 MeV, ro=1.2 fm, and ao=0.687 fm, which yields 
the same barrier height as the Broglia-Winther poten¬ 
tial. The result so obtained is denoted by the solid line 
in the figure. One can see that the one-neutron transfer 
cross sections do not differ much between the two calcu¬ 
lations. In contrast, it is interesting to notice that the 
two-neutron transfer cross sections are considerably af¬ 
fected by the choice of the nuclear potential. Notice that 
we use the same imaginary potential for all the channels. 
Even so, we reach the same conclusion as in Fig. that 
is, the 2n transfer amplitude is more absorbed than the 
In transfer amplitude. In this case, the absorption is 
more likely due to the capture rather than the inelastic 
processes. 



r [fm] 

FIG. 10: (Golor online) The three eigen-barriers obtained 
by diagonalizing the intrinsic Hamiltonian matrix, Eq. fi¬ 
at each internuclear separation, r. 


In order to understand a reasoning for this phe¬ 
nomenon, Fig. shows the eigen-barriers by diagonal¬ 
izing the coupling matrix, Eq. at each position r. 
That is, 

Va{r) = VN{r)^ -l-Aa(r), (11) 

r 

where AQ,(r) is an eigen-value of the coupling matrix. In 
general, the eigen-vectors depend on the position r [SSj. 
However, for the direct two-neutron transfer scheme with 
Voi = Vi 2 and < 52 = 0 , one of the eigen-vectors becomes in¬ 
dependent of r. This special eigen-vector has a structure 
of 1 ^ 1 ) = ^(jOu) — |2n)) with the eigenvalue of — Vo 2 (’’)- 
Notice that the In channel does not contribute to this 
eigen-state. The other two eigen-states are given as a 
liner superposition of all the three channels, jOn), jin), 
and |2n). For a positive value of the ratio /3o2/ao2, that is 
the case in our calculations, the eigenbarrier correspond¬ 
ing to the state {ipi) provides the lowest barrier among 
the three eigenbarriers, and is shown by the dashed line in 
Fig. This implies that the capture occurs more eas¬ 
ily from the 2n channel as compared to the In channel 
due to the lowest eigenbarrier, in which the In channel 
is absent. 

As we will discuss in the next section, the non¬ 
reproduction of the experimental probabilities in the 
vicinity of the barrier shown in Fig. [^may be attributed 
to the fact that we do not include the excited collective 
states in these calculations. 


C. Simultaneous description of fusion and 
multi-neutron transfer 

One of the most important issues in the study of trans¬ 
fer reactions is to investigate whether the same strengths 
for the transfer couplings simultaneously account for fu¬ 
sion and transfer cross sections. Such attempt has been 
successfully made for the 33g_|_90,9i,922j. systems [551 [51] . 
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TABLE VI: The Coulomb deformation parameter, Pc, the nu¬ 
clear deformation parameter, /3jv, the multipolarity and the 
parity, A’^, and the excitation energy E for the collective states 
included in the coupled-channels calculations. The radius pa¬ 
rameter of ro=1.2 fm is used in the coupling potentials. 


Nucleus 

A" 

Pc 

Pn 

E (MeV) 

«Ca 

3" 

0.43 

0.43 

3.737 

^®Zr 

3” 

0.27 

0.305 

1.89 


We do here a similar attempt for the ^°Ca-|-®®Zr reac¬ 
tion, for which fusion cross sections have been measured 
with high precision [ZlIHl. To this end, we include both 
the multi-neutron transfer channels and the collective in¬ 
elastic channels in the coupled-channels equations. To 
be more specific, we take into account the one octupole 
phonon excitation in ^°Ca as well as the octupole phonon 
excitations in ®®Zr up to the three-phonon states. We 
include all the possible mutual excitations. The param¬ 
eters for the collective couplings are given in Table |VI| 
The deformation parameters for the Coulomb couplings 
are estimated with the measured B{E3) values, while we 
slightly increase the nuclear deformation parameter for 
®®Zr in order to better reproduce the experimental data. 
The multi-neutron transfer channels are taken into ac¬ 
count up to three-neutron transfer, with the method of 
Esbensen and Landowne |4]. That is, we treat the collec¬ 
tive excitations and the transfer channels as independent 
degrees of freedom so that the channel wave functions are 
specified as lutrUinei), where ntr (= 0, 1, 2, and 3) indi¬ 
cates the number of transferred neutron while riinei spec¬ 
ifies the inelastic channels. This lead to the total number 
of 32 channels (=2x4x4, where the first 2 is for the 
excitation in ^°Ca, the second 4 is for the excitations in 
®®Zr, and the last 4 is for the transfer channels). For the 
one- and the two-neutron transfers, the parameters for 
the transfer couplings are readjusted to fit the measured 
transfer probabilities with the direct two-neutron transfer 
scheme with the optimum Q-values, that is, Q—0. The 
resultant coupling coefficients are Poi = —71 MeV fm, 
001=1.13 fm, Po 2 = —105 MeV fm, and ao2=0.82 fm. 
The three-neutron transfer channel is included with the 
coupling V 23 =Voi, Vi3=Vb2 and Vb3=0. We employ Qs = 
+4 MeV in order to reproduce the experimental data for 
the three-neutron transfer reaction. We use the nucleus- 
nucleus potential given in Table [Vl 


Figures 11 (a) and 11 (b) show the results of the 
coupled-channels calculations so obtained for the fusion 
cross sections Uf and the fusion barrier distributions Dfus 
for the ^°Ca-|-®®Zr system, respectively. Here, the fusion 
barrier distribution is defined as Df^s = cP{Ea{)/dE‘^ 
[55l[56]. In the figures, the dashed lines denote the re¬ 
sults obtained by including only the inelastic excitation 
in the colliding nuclei, while the solid lines are obtained 
by including in addition the multi-neutron transfer cou- 




FIG. 11: (Color online) (a) The fusion cross sections for the 
reaction. The dashed line shows the results of 
the coupled-channels calculations including only the collective 
excitations in the colliding nuclei, while the solid line shows 
those with both the inelastic and the multi-neutron transfer 
channels, (b) The corresponding fusion barrier distribution. 
The experimental data are taken from Refs. [71[S]. (c) The 
probabilities for the multi-neutron transfer reactions obtained 
with (the solid lines) and without (the dashed lines) taking 
into account the inelastic excitations. The experimental data 
are taken from Ref. [23]. 


plings. One can see that both the fusion cross sections 
and the barrier distributions are well reproduced by in¬ 
cluding the multi-neutron transfer channels, although fu¬ 
sion cross sections are still underestimated at low ener¬ 
gies. In particular, the flatness of fusion barrier distri¬ 
bution is well reproduced by this calculation. We do not 
know the origin for the underestimate of fusion cross sec¬ 
tions at low energies, but a similar tendency is seen also 
in another coupled-channels calculation reported in Ref. 
[46] . This could be due to the absence of other trans¬ 
fer channels, such as proton and alpha particle transfers, 
or the change in the collective couplings for the transfer 
channels. 
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The transfer probabilities are shown in Fig. 11 (c). 
The solid lines show the results of the full coupled- 
channels calculations, while the dashed lines are obtained 
by excluding the inelastic channels, that is, by includ¬ 
ing only the multi-neutron transfer channels. One can 
see that the inelastic excitations do not affect much the 
one- and two-neutron transfer channels, but the three- 
neutron transfer channel is sensitive to the inelastic ex¬ 
citations, especially for D<13 fm. One can also see that 
the present coupled-channels calculations well reproduce 
the experimental data for the transfer probabilities. It is 
worthwhile to mention that the fusion cross sections and 
the transfer probabilities are reproduced simultaneously 
within the single framework. 


IV. SUMMARY 


in order to reproduce the experimental data, ii) the 
higher order dynamics is important at energies around 
the Coulomb barrier and the first order treatment may 
not be sufficient, iii) the absorption has an important ef¬ 
fect on the transfer probabilities at energies around the 
Coulomb barrier, and iv) the absorption plays a more 
important role for the two-neutron transfer channel as 
compared to the one-neutron transfer channel. 

The role of absorption has been confirmed also by using 
the quantal coupled-channels method for the ®°Ni-|-^^®Sn 
system. For the ^°Ca-|-®®Zr system, we have succeeded 
to reproduce simultaneously the fusion cross sections, the 
fusion barrier distribution, and the transfer probabilities 
up to the three-neutron transfer channel, by including 
both the collective excitations and the transfer channels 
in the coupled-channels calculations. 


We have carried out a phenomenological study on 
the heavy-ion multi-neutron transfer reactions using the 
coupled-channels approaches. The aim was to adjust the 
parameters for the transfer couplings using the experi¬ 
mental data at energies far below the Coulomb barrier, 
where the perturbation treatment is applicable, and to 
investigate the dynamics at energies around the barrier, 
where the higher order terms as well as the absorption 
are important. We have applied this strategy to the 
"^°Ca-|-®®Zr and ®°Ni-|-^^®Sn systems, for which the trans¬ 
fer probabilities were recently measured. We have first 
used the semi-classical approximation to a three-channel 
problem with one- and two-neutron transfer channels, 
and have obtained the following conclusions: i) the in¬ 
clusion of the direct coupling between the entrance chan¬ 
nel and the two-neutron transfer channel is necessary 


There have not been many systems for which both the 
fusion and the multi-neutron transfer cross sections are 
available at energies around and far below the Coulomb 
barrier. In that sense, the present phenomenological 
analysis is somewhat limited. It would be interesting 
to test the present approach to future experimental data 
for fusion and multi-neutron transfer reactions. 
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